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The aim of the present work is to investigate the existence of epistatic interactions possibly influencing psychotropic agents' response be- 
tween rs6740584 within Cyclic adenosine monophosphate Response Element Binding (CREB) and rsl2775799 within cAMP response 
element-modulator (CREM) variants in bipolar disorder (BD) and major depressive disorder (MDD). All BD and MDD patients were 
administered with the Young Mania Rating Scale (YMRS) and Hamilton Depression Rating Scale (HAMD) at baseline and at endpoint, 
respectively. A multiple regression model was employed to investigate the existence of possible epistatic interactions between the two 
variants and diverse clinical factors including drug response in affective disorders. No significant epistatic interaction was observed be- 
tween rs6740584 within CREB and rsl2775799 within CREM on both symptom improvement and other clinical factors in affective dis- 
orders. Our preliminary results suggest that no epistatic interaction between rs6740584 within CREB and rsl2775799 within CREM 
should exist on clinical improvement and clinical factors in affective disorders. Psychiatry Investig 2014; 1 1 (2):200-203 
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INTRODUCTION 

An increasing number of studies have recently focused on 
the role of the transcription factor cyclic adenosine monophos- 
phate Response Element Binding (CREB) protein in mood 
disorders such as bipolar disorder (BD). Indeed, CREB plays a 
key role into neuronal signal transduction and it has been im- 
plicated in synaptogenesis, neural connections, and long term 
potentiation. 1 Interestingly, genetic variations within the gene 
coding for CREB (CREB1) have been associated with treat- 
ment outcomes both in patients suffering from major depres- 
sive disorder (MDD) 2 and in patients suffering from BD. 3 In 
addition, CREB is thought to interact with cAMP response el- 
ement-modulator (CREM), a protein that is highly homolo- 
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gous to CREB, with which regulates synaptogenesis and spon- 
taneous neural network activity. 4 Furthermore, CREM has 
been also found to be associated with physiological and devel- 
opmental role within the hypothalamic-pituitary-gonadal axis, 5 
neuronal degeneration and plasticity 6 where are implicated in 
the treatment response of affective disorders. 7 To date, increas- 
ing emphasis has recently been placed upon the investigation 
of possible epistatic interactions between genetic polymor- 
phisms in candidate genes possibly influencing treatment out- 
comes. 8 " 10 Epistasis could be defined as the functional interac- 
tion between genes. It encompasses various events including 
promoter activity control, chromatin remodeling, epigenetic 
control, as well as many other molecular reactions. 11 Because 
these events can impact cell lifecycles and complex traits and 
are orchestrated through genetically driven complex yet flex- 
ible activities, the research for epistatic interactions between 
different candidate genes could represent a significant advan- 
tage in comparison to the investigation of single genes as a way 
to understand the biologic diversity that could influence, for 
instance, drug response. 11 

Accordingly, the aim of the present work is to investigate the 
existence of possible epistatic interactions influencing mood 



200 Copyright © 2014 Korean Neuropsychiatry Association 



AChiesa et al. 



Table 1. Influence of rs6740584 within CREB and rs1 2775799 within CREM on clinical and socio-demographic outcomes in the present 
study 



Clinical and demographic 
characteristics of the sample 



Bipolar 
disorder patients 



Influence of 
rs6740584 



Influence of 
rs!2775799 



Influence of the epistatic 

interaction between 
rs6740584 and rs!2775799 





^1M — 1 JZ J 


P 


p-value 


B 


p-value 


P 


p-value 


Gender 
















Males 


87 (66%) 


0.42 


0.07 


0.46 


0.03 


-0.54 


0.05 


Age 


36.35±11.60 


0.14 


0.55 


0.06 


0.76 


-0.24 


0.38 


YMRS baseline 


33.27±9.09 


-0.15 


0.51 


-0.03 


0.86 


0.02 


0.94 


YMRS discharge 


19.80±5.25 


-0.08 


0.72 


0.04 


0.83 


-0.03 


0.90 


YMRS % improvement 


37. 14± 19.26 


0.08 


0.72 


0.13 


0.49 


-0.13 


0.63 


Response 
















Yes 


33 (25%) 


0.24 


0.30 


0.27 


0.17 


-0.27 


0.33 


Remission 
















Yes 


10 (8%) 


0.08 


0.72 


0.07 


0.73 


-0.07 


0.78 


Clinical subtypes 
















BD without PF 


73 (55%) 


-0.40 


0.08 


-0.34 


0.07 


0.59 


0.04 


BDwithPF 


57 (43%) 














Age at onset 


26.58±10.19 


0.39 


0.11 


0.29 


0.14 


-0.54 


0.07 


Familiar history of psychiatric disorders 
















Yes 


42 (32%) 


0.19 


0.52 


-0.18 


0.46 


0.23 


0.50 


Suicide attempts 
















Yes 


22 (17%) 


-0.19 


0.50 


-0.10 


0.62 


0.26 


0.41 


Duration of admission (days) 


33.66±21.07 


0.02 


0.93 


0.14 


0.50 


-0.11 


0.71 


Drug 
















Lithium 


41 (31%) 


0.17 


0.54 


0.18 


0.41 


-0.26 


0.43 


Valproate 


56 (42%) 














Other 


5 (4%) 














Concomitant anxiolytics 
















Alprazolam 


5 (4%) 


0.35 


0.19 


0.25 


0.26 


-0.41 


0.21 



Lorazepam 
Clonazepam 
Buspiron 
None 



82 (62%) 
3 (3%) 
7 (5%) 

35 (26%) 



CREB: cyclic adenosine monophosphate Response Element Binding, CREM: cAMP response element-modulator, YMRS: Young Mania Rat- 
ing Scale, BD: bipolar disorder, PF: psychotic features 



stabilizers and antidepressants response between genetic vari- 
ants within CREB and CREM that we individually investigated 
in a previous paper 12 in BD and MDD patients treated with 
different mood stabilizers and antidepressants. 

METHODS 

The subjects was described in detail in a previous study. 12 
Briefly, 132 and 145 in-patients suffering from BD and MDD 
according to DSM-IV criteria 13 were consecutively collected at 
the Department of Psychiatry of the Catholic University of 
Korea College of Medicine, Seoul, Korea. There was no partic- 
ular restriction with regard to the specific mood stabilizers em- 



ployed. However, patients were excluded if they had current se- 
vere or unstable medical and neurological conditions, current 
treatment with a long-acting antipsychotic, concomitant alco- 
hol and substance abuse disorders and if they were not of Ko- 
rean ethnicity. All patients admitted to the hospital were as- 
sessed for the severity of illness at baseline and at discharge 
by means of the Young Mania Rating Scale (YMRS) and Ham- 
ilton Depression Rating Scale (HAMD). 14,15 In accordance with 
previous studies, response was a priori defined as a >50% 
symptoms' reduction from baseline to discharge. 16 Remission 
was denned as a YMRS score <12 and HAMD score <7 at dis- 
charge, respectively. 17 As an example, Table 1 represent the 
case of bipolar patients. The case of MDD is not presented due 
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to space limitation (data available on request). The study pro- 
tocol was approved by the institutional review board (approval 
number HC10TISI0031). 

The main outcome measure of the present study was the ex- 
istence of possible epistatic interactions between rs6740584 
within CREB [position=208429351(8922), alleles=C/T, loca- 
tion=intron)] and rsl2775799 within CREM {[position= 
35488264 (61501), alleles=C/T, location=intron]} and YMRS 
and HAMD improvement scores. Secondary outcome mea- 
sures included the investigation of epistatic interactions be- 
tween the same two SNPs and other clinical and socio-de- 
mographical variables included in the present study. The two 
SNPs mentioned above were selected in accordance with the 
following criteria: 1) among SNPs available in our dataset, 12 2) 
we chose 1 SNP for each gene, 3) that tagged the largest num- 
ber of SNPs within the same gene, so as to cover a portion of 
the gene as high as possible (data from http://hapmap.ncbi. 
nlm.nih.gov/). For CREB, only one SNP was available. 12 For 
CREM, we choose rsl2775799 because it was the SNP that 
allowed us to cover the larger part of the gene (approximately 
14%). 

Statistical analyses were performed using 'Statistica pack- 
age. 18 A multiple regression model was employed to investigate 
the existence of possible epistatic interactions between the two 
genotypes and clinical/socio-demographical variables includ- 
ed in the present study. Clinical improvement was calculated 
based on changes from baseline to the end of treatment. 

All p-values were 2-tailed and statistical significance was 
conservatively set at the 0.003 level (approximately corre- 
sponding to the 14 variables outlined in Table 1). With these 
parameters we had a sufficient power (0.80) to detect a medi- 
um-large (d=0.3) effect size for patients carrying the TT rs 
6740584 genotype as compared with those carrying the CT 
genotype. 19 Such effect size corresponded to the possibility of 
detecting final differences on YMRS as high as 2 points. 

RESULTS 

No significant epistatic interaction between rs6740584 with- 
in CREB and rsl2775799 within CREM on YMRS improve- 
ment measured by the change in YMRS total scores from 
baseline (admission) to the endpoint (discharge) was observed 
((3=-0.13, p=0.63). In the analysis of secondary outcome mea- 
sures (total score on YMRS at baseline and discharge, re- 
sponse, remission, subtypes of BD, age, gender, age at onset, 
family history of psychiatric disorders, duration of admission, 
suicide attempt, mood stabilizers, and concomitant anxiolyt- 
ics), no significant epistatic interactions between rs6740584 
within CREB and rsl2775799 within CREM were observed as 
well (all p-values >0.003). Similar trends were also observed in 



HAMD -related outcomes and clinical factors in patients with 
MDD. A multiple regression results are presented in Table 1 
(BD data only and MDD data available on request). 

DISCUSSION 

The present work was aimed at exploring the existence of 
possible epistatic interactions between rs6740584 within 
CREB and rsl2775799 within CREM on symptom improve- 
ment scores and on other clinical variables in a sample of BD 
and MDD patients. The results of our study point to the lack 
of significant epistatic interactions between the two SNPs and 
clinical improvement, suggesting a lack of association be- 
tween these genetic variants and clinical improvement in BD 
and MDD patients. Note, however, that our finding could also 
be due to the moderately small sample size of our sample that 
could not allow us to detect subtle differences that are usually 
associated with single genes or gene-gene interactions in com- 
plex disorders. 20 

It is also worth mentioning that in our earlier paper focus- 
ing in single SNPs within CREB and CREM no significant as- 
sociation with YMRS scores was observed as well. 12 However, 
there is some evidence so far to suggest that epistatic interac- 
tions can be observed even when single SNPs individually 
considered do not lead to positive findings. 21 Such a possibility, 
in turn, raises significant computational and theoretical diffi- 
culties when considering all possible gene-gene interactions 
not to mention other forms of genetic expression controls. 22 
Note also that no further significant epistatic interactions on 
any of the secondary outcome measures have been observed. 
These findings are in line with earlier findings observed in our 
original study. 12 Moreover, they extend those findings by sug- 
gesting that no epistatic interaction exists between the two 
SNPs. Not surprisingly, no interaction effects between CREB 
and CREM on phenotypes and treatment outcomes were also 
seen in patients with MDD in our further analysis. Hence, we 
may assume that these two variants do not have epistatic in- 
teraction at least in our patients with affective disorders. 

Before firm conclusions are drawn, several limitations af- 
fecting the present study should be considered. 23 First of all, 
the lack of associations observed in the present study could 
be simply due to the lack of statistical power. This, in turn, 
could obscure small effects exerted by single SNPs. This issue 
is particularly concerning if one considers that even among 
pharmacogenetic studies with a large sample size, results tend 
to remain conflicting. 24 A further concern is related to the use 
of mood stabilizers with different mechanisms of action that 
do not allow to reach definitive conclusions with regard to the 
influence of the SNPs on single mood stabilizers. 25 However, 
our decision to include patients treated with different drugs 
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could have the advantage of being closer to "real world" clini- 
cal practice. In addition, the duration of hospitalization in the 
present study could be considered as insufficient to ascertain 
a lack of response and remission. However, this time frame is 
consistent with common clinical practice. 26 Furthermore the 
different duration of hospitalization could raise concerns about 
the fact that clinical improvement could vary as a function of 
time rather than of genetic variants. However, we checked for 
such a possibility, finding no significant influence. Finally, a 
further limitation of this study could be related to the incom- 
plete coverage of genes due to our tagging approach. 

Epistasis between different genes still represents a particu- 
larly relevant analysis to be performed in genetic association 
studies till today, being helped by the knowledge of the prin- 
cipal ways by which the so-called reactome displays its activi- 
ties. 9 ' 10 The inclusion of the analysis of the epistasis between 
genes along with the analysis of the impact of single genotypes 
and haplotypes should be considered as a standard in the next 
association studies. 9,10 When epistasis findings are widely 
available and collected across research teams over the world, 
such data will also substantially contribute to the understand- 
ing of the role of gene-gene interaction in affective disorders 
and establishment of personalized-medicine for difficult-to- 
treat patients. However, our findings preliminary suggest that 
no epistatic interactions exist between CREB and CREM and 
clinical outcomes in BD and MDD patients treated with dif- 
ferent antidepressants. Thus, further researches may investi- 
gate whether our results hold in a larger sample of BD and 
MDD patients treated with different antidepressants or wheth- 
er the investigation of different SNPs and/ or different gene- 
gene interactions within other genes involved in the transcrip- 
tion pathways might lead to different outcomes. 
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